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ABSTRACT 

We analyze deep X-ray, optical and mid-infrared Spitzer observations of the CDF-N/GOODS-N 
region to study a sample of 13 submillimeter-detected galaxies with spectroscopic redshifts (median 
2 = 2.2). These galaxies are among the most active and massive at this epoch. We find evidence for 
a power-law correlation between the estimated stellar and X-ray luminosity, implying that masses of 
the black holes may be related to the stellar masses of their host galaxies. We derive the rest-frame 
UV-near-infrared spectral energy distributions for these galaxies, believed to be young spheroids, and 
fit them with model templates. Using the rest-frame near-infrared luminosities, which are relatively 
insensitive to uncertainties in stellar ages and reddening in these young dusty galaxies, and theoretical 
mass-to-light ratios, we can estimate their stellar masses. Although the submillimeter emission implies 
that these galaxies are undergoing an epoch of intense star-formation, the Spitzer data reveal a massive 
stellar population already in place. We find that our submillimeter galaxies have a median stellar mass 
of ~ 10"'^^ M0, which is roughly ten times more massive than typical UV-selected star-forming systems 
at similar redshifts. These stellar masses are then compared to previously published black hole mass 
estimates derived from the X-ray luminosities under the assumption of Eddington-limit accretion. 
We find that the black hole masses for our high-redshift sample are approximately 1-2 orders of 
magnitude smaller than galaxies of comparable stellar mass in the local Universe. Although our 
estimates of black hole masses will increase if the accretion is sub-Eddington, and our stellar masses 
will decrease if we assume a much younger stellar population or a different initial mass function, we 
find that only through a combination of effects is it possible to shift the high redshift galaxies such 
that they lie on the local relation. This suggests that the black holes need to grow substantially 
between z = 2.2 and the present-day, with much of the black hole growth occurring after the current 
obscured, far-infrared luminous phase of activity which is likely associated with the formation of 
the spheroid. This interpretation supports a scenario where submillimeter galaxies pass through a 
subsequent accretion-dominated phase, where they would appear as optically bright quasars. 
Subject headings: galaxies: evolution - galaxies: formation - galaxies: starburst 



1. INTRODUCTION 

Studies of galaxies in the local Universe have demon- 
strated an apparent relation between the stellar mass 
of a spheroid and that of the super-massive black hole 
(SMBH) which hes at its heart. The M*-Mbh relation 
for spheroids, with the mass of the SMBH typically ~ 
0.1% of the spheroid r nass ( e.g. iKormeiKbL-fcJ^i ahstong 
119951 rMagorri ~et al.lll99at iGebhardt et al.ll200(t . s~ 
gests that the growth of the black hole and the surround- 
ing galaxy are related. This provides an important clue 
to understanding galaxy formation, as SMBHs have been 
suggested as a highly efficient mechani sm for regulat- 
ing star formation in galaxies (e.g . Caval iere fc Vittorinil 
I200I iGraiiato et al.ll200l: IDi Matteo et al.ll2005fl . Such 
feedback is an essential part of current galaxy forma- 
tion models and provides a physically-motivated expla- 
nation for the form of the galaxy lu minosity function 
(|Benson et al.lf 2003: Ba ugh et al.]l20n,5D . However, obser- 
vational evidence for feedback, in particular driven by ac- 

Electronic address: borys@caltech.edu 

1 Caltech, 1200 E California Blvd., Pasadena, CA 91125 

^ Institute for Computational Cosmology, Durham University, 
South Road, Durham DHl 3LE, UK 

^ Institute of Astronomy, University of Cambridge, Madingley 
Road, Cambridge CB3 OHA, UK 

Astronomy Technology Centre, Royal Observatory, Blackford 
Hill, Edinburgh EH9 3HJ, UK 

^ Institute for Astronomy, University of Edinburgh, Blackford 
Hill, Edinburgh EH9 3HJ, UK 



tivity from SMBHs, is sparse, especially for the youngest 
and most massive g alaxies - for which it is expected to 
be most important iBenson et alJl2(¥)^ . 



3nson 

Modeling AGN-driven feedback is more difficult than 
the traditional star-formation feedback recipes. In part 
this is because we still lack a clear model for the early 
growth of the SMBH, where BH-BH mergers and super- 
Eddington accretion may both play a part. There is thus 
an urgent need for observations of the earliest phases of 
the evolution of spheroids and SMBHs to guide develop- 
ment of theoretical models. 
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sive galaxies. The intensity of their 
resulting high metallicity, along with their 
na mical masses, high gas fractions and strong cluster- 
ing flvison et al."2002';' Greve et alJ2005tlSwinbank et alJ 
1200 4; Blain et al. 2004) are all suggestive of a close link 
to the formation phase of the most massive spheroids 
llSmail et alJl200l . If these are indeed young spheroids, 
what can we learn about their black holes? 

Deep X-ray observations provide a sensitive probe of 
AGN activity in submillimeter galaxies (SMG) and can 
be used to place SMBH mass constraints, even in the 
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presence of heavy gas and dust obscuration. Exploit- 
ing the unique sensitivity of the 2-Ms exposure of the 
Chan d ra Deep Fie ld Nort h (CD F-N; [Alexander et all 
12003ft . [Alexander e t al. (20 05albD study a small sample 
of spectra scopically identified SMGs and find that most 
host weak AGN with modest absorption-corrected X- 
ray luminosities. This suggest two things: 1) the AGNs 
do not typically dominate the bolometric emission from 
SMGs, and 2) the masses of the SMBHs are likely to be 
< IO^Mq (Eddington-limited SMBH masses are typically 

IO'^Mq). The large AGN fraction in the SMG popula- 
tion (>40% for sources with /sso^m ~ 4 mJy) further in- 
dicates that the SMBHs are growing almost continuously 
throughout these vigorous, but obscured, star-formation 
episodes. The stark contrast between the AGN fractions 
in the SMG and coeval U V-bright galaxy populations 
(~ 8%: lSteidel et al.ll|2004D ) suggest that this joint black 
hole-stellar growth activity is most strongly connected 
with the SMG population (Alexander et al. 2005a). 

In this paper we compare estimates for the stellar and 
SMBH masses for a sample of 13 high-redshift proto- 
spheroids (identified with far-infrared luminous galaxies 
detected in the submillimeter and radio wavebands) in 
the GOODS-N region. These all have precise spectro- 
scopically identified redshifts and have hard X-ray emis- 
sion dominated by an AGN. The stellar mass estimates 
are derived from optical, near-infrared and most impor- 
tantly the mid-infrared observations of this region by the 
Syitzer S pace Telescope tak en as part of the GOODS sur- 
vey (Dickinson et al. 2005), whi le the S MBH mass esti- 
mates arc taken from Alexan der et alJ f2QQ5a, b) based 
on the 2-Ms Chandra observations. Together, these data 
allow us to investigate the origin of the M^^-Mbh rela- 
tion for young spheroids in the early Universe. Through- 
out, we adopt a cosmology with Hp ~ 65kms~^ Mpc~^, 
^^M — 1/3, and Ha = 2/3. Unless otherwise stated, 
masses and luminosities are given in solar units. 

2. SAMPLE SELECTION AND DATA REDUCTION 

For l ocal galaxies, spatially resolved kinematics (e.g. 
M^p;9rrian et al. 1998) or reverberati on mapping (e.g. 
Wandel et aP Il999l: iKaspi et all 12000(1 are useful tech- 
niques for weighing the SMBHs at the center of galaxies 
and relating these to the stellar masses of their spheroids. 
However, neither of these techniques can currently be ap- 
plied to SMGs at high redshifts and so by necessity we 
have to adopt much cruder methods to estimate both 
the black hole and spheroid masses. In this study we use 
the emission arising from accretion onto the SMBH as 
a tracer of the black hole mass and the rest-frame near- 
infrared luminosity of galaxy as an indicator of its stellar 
mass. 

Unfortunately, the rest-frame UV and optical spectra 
of the SMGs provides an ambig uous tracer of the nu- 
clear activity in these galaxies ('Ghao man et al.l 120051: 
l&winbank et al. 2004). The difhculty of merely identi- 
fying an AGN from these observations, let alone reliably 
estimating their luminosities is likely caused by heavy 
obscuration towards the SMBH or faint broad lines from 
intrinsically weak AGN being lost against the continuum 
of the galaxy. Conversely the weakness of the AGN and 
their obscuration makes the use of the rest-frame near- 
infrared luminosity a much more reliable measure of the 
stellar mass of the galaxies. 



In contrast, hard X-ray emission from the accre- 
tion disk emerges relatively unattenuated from all but 
the most obscured regions, hence the X-ray emission 
from SMGs may provide the most reliable probe of 
their SMBH ma sses. However, as demonstrated by 
[Ale xa nder et all l|2005a(l . only the Chandra 2-Ms ob- 
servation of the CDF-N region ((Alexander et al.l 12003(1 
reaches the necessary depth to provide constraints on 
the X-ray luminosities of typical SMGs. For this reason 
our analysis focuses on the sam ple of SMGs with reli- 
able spectroscopic redshifts from IChapman et alJ ((20031 
12005ft which lie within this field. We stress that pre- 
cise redshifts are essential for our analysis. Not only are 
the mass estimates dependent on a well-constrained lu- 
minosity distance, but more critically there are strong 
degeneracies in redshift and age in the fitting of their 
spectral energy distributions (SED) described in §3. 

2.1. X-ray constraints on SMGs 

Using the sample of spectroscopically identified, radio- 
detect ed SMGs in the CDF-N from IChapman et alJ 
l(2flO;^ , obscuration-corrected X-ray luminosities were 
derived for 15 galaxies and SMBH masses esti- 
mated assuming Eddington accretion by Alexander et alJ 
|2005albl ) . We discard two of these galaxies from our 
analysis: J123716.0H-620323.3 falls outside the mid- 
infrared coverage, and the optical/mid-infrared photome- 
try of J123632. 61+620800.1 is contaminated by a diffrac- 
tion spike from a bright {Kg ~ 14) object only 11" away. 
These conditions prevent us from estimating the stellar 
masses needed in §3.2, hence we concentrate on 13 SMGs 
with optical, near-infrared and mid-infrared photome- 
try or limits and reliable spectroscopic redshifts. These 
galaxies range in redshift from z = 0.5-2.9, with a me- 
dian comparable to that of the entire SMG population 
{z = 2.2). 

The X-ray spectral analysis in lAlexander et aTl l(2005bfl 
shows that the bulk of the hard X-ray emission from 
the SMGs in our sample is likely to be due to accre- 
tion onto the SMBH, with star formation contributing 
only at the softest energies, ^ 2keV in the rest-frame. 
The luminosity from the X-ray waveband is only a few 
percent of the bolometric luminosity from the SMBH, 
and [Ale xander et al. ( 2005b.) adopted a bolometric cor- 
rection (BCx) of 6^4^% consistent with other studies. 
They also employed Eddington arguments to show that 
log(MBH) = log(iboi/??) — 38.06, where the bolomet- 
ric luminosity is in units of ergs^^. The factor rj is 
unity for Eddington-limited accretion, and < 1 for 
sub-Eddington. As our bench-mark we adopt Eddington- 
limited accretion and use this to investigate the relation 
between black hole and stellar masses for our sample. 

2.2. Spitzer observations 

In addition to the X-ray constraints on the luminosi- 
ties and hence masses of the SMBHs, our analysis also 
requires estimates of the stellar masses of the host galax- 
ies. A particularly simple and relatively robust way of 
deriving the mass of a stellar population is to use the rest- 
frame near-infrared luminosity and a theoretical mass 
to light ratio. For all but the youngest stellar popu- 
lations, the rest-frame near- infrared emission is domi- 
nated by light from main sequence stars and hence, for 
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example, the rest-frame if-band luminosity changes by 
only a factor of two bet ween stellar populatio ns with ages 
of 100 Myr and 1 Gyr ijLeitherer et alJll999() . The rest- 
frame y-band luminosities of the same populations on 
the other hand vary by a factor of ten. In addition to 
its relative insensitivity to age variations, the rest- frame 
-fsT-band is over ten times less sensitive to reddening than 
V-band, which is important to consider for these dusty 
SMGs. 

At the redshift of our targets, the rest frame i^T-band 
is shifted towards the mid-infrared and hence Spitzer ob- 
servations are required to estimate the stellar mass. For 
this reason, we obtained the IRAC images of GOODS- 
N from the Spitzer archive (V.l 20041027 and V.l 
20050505 enhanced data products). We used Sextrac- 
TOR (Bertin & Arnouts 199^ to measure the photome- 
try and uncertainties in these for each of the galaxies on 
the IRAC imaging. By constructing curves of growth for 
several isolated objects, we determined that 6" diame- 
ter apertures were necessary to measure the total IRAC 
flux from our targets. Since source confusion is an issue 
for these extremely deep exposures, we used 4" diameter 
apertures and then applied a small aperture correction, 
determined to be 0.22, 0.26, 0.33, and 0.30 mag for the 
3.6, 4.5, 5.8, and 8.0 /im bands respectively. Even so, we 
were unable to adequately constrain the 3.6 /im flux for 
three sources due to confusion noise from nearby sources. 
We provide the IRAC photometry for all 13 sources in 
Table ni At 5.8/im, which corresponds to a rest-frame 
wavelength of ~ 2/im, the ten high-redshift SMGs at 
z > 1.8 have a bright mean magnitude of 20.1 ±0.8, with 
the entire sample spanning under an order of magnitude 
in flux. This is in contrast with the i?-band fluxes, which 
cover well over an order of magnitude. The modest dis- 
persion in the mid-infrared fluxes of the SMGs already 
suggests that we are dealing with a relatively homoge- 
neous population of luminous, high redshift galaxies. 

2.3. Spectral Energy Distributions 

To improve the constraints on the form of the SEDs in 
SMGs we also make use of rest-frame UV/optical pho- 
tometry to extend the wavelength coverage and so esti- 
mate the likely reddening and ages of the stellar popu- 
lations. For this we use the UBVIRz' photometry pre- 
sented in iCaoak et al.1 l)2004f ). These data were supple- 
mented with J- and i^g-band near-i nfrared imaging frorn 
the WIRC2 cani e ra at Palomar (see lBundv et all l^2005^ 
and ISmail et"an (|2004[) for a description of these data) . 
Combined, these ground-based data cover more bands 
than the GOODS-N Hubble Syace Telescope (HST) ACS 
imagmg of this field ijGiavalisco et al.ll2004D . Although 
the extra bands provided by the ground based imaging is 
essential for the SED template fitting (see ii3.1|l , the high- 
resolution imaging available from HST is valuable for 
studyi ng the light profiles of these galaxies. iBorvs et al.l 
l)2006t) show that the rest-frame F-band light distribu- 
tions of a sample of 27 SMGs (5 of which are in the 
sample presented here) follow a r^^^ power law on aver- 
age, consistent with the spheroidal-nature of these sys- 
tems. A summary of the optical and near-IR photometry 
is presented in Table |21 We adopt 3" and 4" diame- 
ter apertures to measure the fluxes from the optical and 
near-infrared imaging as these provide optimal signal-to- 
noise and a near-total estimate of the magnitudes for the 



SMGs. 

We plot the combined optical, near- and mid-infrared 
photometry for the 13 SMGs in Fig. [T] The SEDs for 
these galaxies show three clear features: a sharp decline 
at rest-frame wavelength of < 912A, resulting from the 
Lyman— a limit within the galaxy and the foreground 
intergalactic medium; a break at around 4000A aris- 
ing either from the spectral breaks found in young or 
evolved stellar populations at this approximate wave- 
length and finally a peak in the emission at ^ 1.6/im 
associated with the mi nimum in t he opacity ( Johnl 
[T988llSimns~fc Eisen harddll999D and hence the maxi- 
mum emission from stellar atmospheres. 

Looking at the SEDs in Fig. we conclude that there 
is a modest range in the rest- frame UV/optical/near- 
infrared properties of radio-detected SMGs. However, 
focusing on the ten higher-redshift SMGs, at z > 1.8, of 
those with full coverage of their rest-frame UV and op- 
tical SEDs, 80% show a detectable discontinuity in their 
SEDs around 4000A in the rest-frame (Fig. QJ. We in- 
terpret this break as arising from the Balmer limit in 
stellar populations with ages of ~ lOO's Myrs. This in- 
terpretation is consistent with the likely nature of SMGs 
as luminous far-infrared sources seen during a phase of 
intense star formation. We note that age-related differen- 
tial reddening may also account for the relative strength 
of this feature in an even younger on-going starburst 
iPoggianti fc Wdl200(il) . 

3. ANALYSIS AND RESULTS 

3.1. Estimating Stellar Masses 

Due to the complex mix of extinction expected within 
SMGs ijChapman et al.ll2004() . which may influence dif- 
ferent age stellar populations to different degrees, we 
have adopted a very simple approach in determining 
our stellar masses. We first estimate the rest-frame K- 
band luminosities of the galaxies, interpolated if neces- 
sary from a fit to the SED, and then combine these with 
estimates of the mass-to-light ratios in the iiT-band for a 
range of plausible ages for the dominant stellar popula- 
tion. 

To derive the rest-frame -f C-band luminosities we ex- 
ploit the HYPER-z package ijBolzonella et all \200(t) to 
fit physically meaningful spectral templates to the spec- 
tral energy distributions of each galaxy and then esti- 
mate their A'-band luminosities by interpolation of the 
best-fit templates. This approach has the added ad- 
vantage that the properties of the best-fit spectral tem- 
plates may provide some insights into the properties of 
the stellar populations within these galaxies. We restrict 
the code to fit the SEDs at the spectroscopically con- 
firmed redshift, and allow it to choose between two model 
SED templates: one having an instantaneous burst of 
star-formation, or one with a constant star formation 
(CSF) history. These represent the two extremes of the 
likely star formati on histories for these young galaxies 
(jSmail et alJl200^ and provide an adequate variety to 
guarantee a reliable estimate of the K-hand luminosity 
for those galaxies whose SEDs are dominated by star- 
light. In our fits, the reddening (Ay) was free to vary, 
and the ages at the galaxy's redshift are constrained to 
be physical within our adopted cosmology. These provide 
additional degrees of freedom for the fitting and may give 
some indication of the properties of the composite stellar 
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TABLE 1 
Spitzer IRAC photometry 



ID 


3.6/im 


4.5/im 


5.8^m 


8.0/im 


J123636.75+621156.1 


20.61 + 0.08 


20.96 + 0.08 


21.18 + 0.22 


21.44 + 0.14 


J123721. 87+621035. 3 


19.35 ± 0.04 


19.66 + 0.04 


19.96 + 0.11 


20.12 + 0.05 


J123629.13+621045.8 


19.06 + 0.04 


19.13 + 0.04 


19.37 + 0.10 


19.59 + 0.06 


J123555.14+620901.7 


19.91 + 0.10 


19.51 +0.08 


19.31 + 0.17 


19.29 + 0.09 


J123711. 98+621325. 7 


21.56 + 0.10 


21.33 + 0.08 


21.08 + 0.15 


21.52 + 0.08 


J123635.59+621424.1 


19.45 + 0.04 


18.98 + 0.03 


18.46 + 0.05 


18.03 + 0.02 


J123549.44+621536.8''''' 




19.78 + 0.30 


19.61 + 0.30 


19.63 + 0.30 


J123606.72+621550.7 


21.52 + 0.17 


21.18 + 0.11 


20.85 + 0.25 


20.68 + 0.11 


J123622.65+621629.7 


20.74 + 0.08 


20.38 + 0.05 


20.19 + 0.11 


20.69 + 0.06 


J123707.21+621408.1 


20.78 + 0.10 


20.44 + 0.08 


20.18 + 0.14 


20.55 + 0.05 


J123606.85+621021.4 


20.37 + 0.08 


20.25 + 0.06 


20.19 + 0.15 


20.61 + 0.10 


J123616.15+621513.7i^ 




20.76 + 0.08 


20.36 + 0.15 


20.05 + 0.07 


J123712.05+621212.3^ 




21.06 + 0.07 


20.78 + 0.13 


20.83 + 0.05 



Note. — All magnitudes are in the AB system. Values represent 6" diameter apertures, derived by applying a correction to the measured 
4" aperture magnitudes. 

^This sources lies on the edge of the image, although is clearly detected. In these cases we set the error to 0.3 mag to reflect the uncertainty 
in our aperture fluxes. 

''The 3.6/xm images for these sources were confused, and sextractor could not converge on a solution. Since the other photometry is 
sufficient to constrain the shape of the SED, we choose not to use the data from this band. 



TABLE 2 

Ground-based optical and near-IR photometry 





ID 




U 








B 








V 








R 








/ 








z' 






J 






Ks 


a 




J123636, 


.75+621156, 


1 


24.28 + 


0, 


,08 


23, 


,43 + 


0, 


,05 


23, 


,09 + 


0, 


,03 


22, 


.39 + 


0, 


,02 


21, 


,84 + 


0, 


,02 


21, 


,61 + 


0, 


,03 


21.12 + 


0, 


,04 


20.58 + 


0, 


,04 


J123721, 


.87+621035, 


,3 


24.19 + 


0, 


,08 


24, 


,00 + 


0, 


,07 


23, 


,52 + 


0, 


,04 


22, 


,81 + 


0, 


,02 


21, 


,84 + 


0, 


,03 


21, 


,35 + 


0, 


,02 








19.84 + 


0, 


,02 


J123629, 


,13+621045, 


,8 


26.11 + 


0, 


,25 


25, 


,70 + 


0, 


,15 


24, 


,98 + 


0, 


,09 


24, 


.05 + 


0, 


,04 


22, 


,89 + 


0, 


,04 


22, 


,36 + 


0, 


,04 


21.39 + 


0, 


,05 


19.93 + 


0, 


,02 


J123555, 


,14+620901, 


,7 


24.38 + 


0, 


,09 


24, 


,24 + 


0, 


,08 


24, 


,16 + 


0, 


,05 


23, 


.80 + 


0, 


,04 


23, 


,45 + 


0, 


,06 


23, 


,16 + 


0, 


,06 


21.99 + 


0, 


,09 


20.92 + 


0, 


,05 


J123711, 


,98+621325, 


,7 


26.12 + 


0, 


,26 


25, 


,75 + 


0, 


,14 


25, 


,71 + 


0, 


,16 


25, 


,66 + 


0, 


13 


24, 


,99 + 


0, 


,18 


25, 


,06 + 


0, 


,24 


> 24.0 






> 23.1 






J123635, 


,59+621424, 


1 


23.94 + 


0, 


,07 


24, 


,06 + 


0, 


,07 


24, 


,00 + 


0, 


,05 


23, 


.79 + 


0, 


,04 


23, 


,36 + 


0, 


,06 


23, 


,11 + 


0, 


,06 


21.99 + 


0, 


,09 


20.74 + 


0, 


,04 


J123549, 


,44+621536, 


,8 


24.61 + 


0, 


10 


23, 


,93 + 


0, 


,07 


23, 


,60 + 


0, 


,04 


23, 


,27 + 


0, 


,03 


23, 


,01 + 


0, 


,04 


22, 


,84 + 


0, 


,05 


21.77 + 


0, 


,07 


21.49 + 


0, 


18=^ 


J123606, 


,72+621550, 


,7 


23.75 + 


0, 


,06 


23, 


,32 + 


0, 


,06 


23, 


,29 + 


0, 


,03 


23, 


.26 + 


0, 


,03 


23, 


,09 + 


0, 


,05 


23, 


,06 + 


0, 


,06 


> 24.0 






> 22.8^ 






J123622, 


,65+621629, 


,7 


25.62 + 


0, 


17 


25, 


,17 + 


0, 


,14 


25, 


,09 + 


0, 


,11 


24, 


.85 + 


0, 


,08 


24, 


,45 + 


0, 


,13 


24, 


,21 + 


0, 


,14 


> 24.0 






22.25 + 


0, 


15 


J123707, 


,21+621408, 


1 


> 27.8 






26, 


,77 + 


0, 


,27 


26, 


,36 + 


0, 


,29 


25, 


,63 + 


0, 


14 


25, 


,51 + 


0, 


,30 


24, 


,69 + 


0, 


,19 


> 24.0 






21.86 + 


0, 


11 


J123606, 


,85+621021, 


,4 


26.17 + 


0, 


,26 


25, 


,25 + 


0, 


,14 


25, 


,15 + 


0, 


,10 


24, 


,41 + 


0, 


,06 


24, 


,21 + 


0, 


,10 


23, 


,93 + 


0, 


,10 


22.20 + 


0, 


10 


20.98 + 


0, 


,05 


J123616, 


,15+621513, 


,7 


> 27.8 






26, 


,59 + 


0, 


,24 


26, 


,13 + 


0, 


,22 


25, 


,36 + 


0, 


11 


25, 


,03 + 


0, 


,18 


25, 


,00 + 


0, 


,23 


> 24.0 






22.71 + 


0, 


,23 


J123712, 


,05+621212, 


,3 


> 27.8 






> 27.6 






> 27.6 






> 27.4 






> 26.4 






> 26.2 






> 24.0 






23.02 + 


0, 


,30 



Note. — All magnitudes arc in the AB system. For undetected sources, we quote 2.5ct upper limits. The J and Kg magnitudes are measured in 
4 apertures, while all other bands use 3 . 

^Two sources that fall outside our Ks coverage are in the HK' catalog provided bv lCapak et al so wc use those values here. 



population in these systems. 

For each galaxy, hyper-Z returns the absolute magni- 
tude in the rest-frame if -band, as well as the age, redden- 
ing and star formation history of the best-fit template. 
Although the M-k values are not corrected for it, the red- 
dening (whic h is ^ 0.1 Ay based on the extinction law 
presented in iCalzetti et alJ 1)200(11) ) has a negligible ef- 
fect on our stellar mass estimates, particularly compared 
with other uncertainties in our analysis. The typical fit 
to a SED in our sample produces a reasonable x^, with 
only three minor exceptions: Both J123635. 59+621424.1 
and J123606. 72+621550. 7 are spectroscopically identi- 
fied AGN, and demonstrate an obvious excess at the 
reddest wavelengths. This would be expected, since 
the input templates do not account for the contri- 
bution in the mid-IR of warm dust from an AGN. 
J123549. 44+621536. 8 also has a moderate near-IR ex- 
cess, but the significance is muted somewhat by its 
proximity to the edge of the IRAC map, where pho- 



tometry was more difficult to estimate. The starburst 
J123606.85+621021.4 has an excess in the UV, but this 
is caused by a contribution from an unrelated blue fore- 
ground {z = 0.4) galaxy only 2" away. In all these cases, 
we feel the best-fit SED template is still a good measure 
of the rest-frame if-band luminosity and provides a good 
estimate of the underlying stellar mass. As we show later, 
changes to the implied stellar masses of factors of 2 (due 
to AGN contamination at restframe if-band), would not 
significantly alter the conclusions of this work. 

To convert these estimates for Mif into stellar masses 
we need to determine the light-to-mass ratio, L^f/M for 
the dominant stellar populations in these galaxies. To 
do this we tur n to the Starburst99 stellar population 
model (Leith crer et al.l[l999,) , which provides estimates 
of for our two extreme scenarios - either a burst or 
constant star formation hist ory. We adop t the the ini- 
tial mass function (IMF) of IMiller fc Sca le (1979), this 
is also the same as that used for the hyper-z templates. 
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Fig. 1. — The spectral energy distributions of the 13 SMGs with black hole mass estimates from I Alexander et alJ I2n05bl) . ordered in 
increasing redshift. The top line shows the SEDs of the lower redshift SMGs, while the remaining 10 are all z > 1.8. The best-fit template 
from HYPER-z is overlaid on the measured photometry (in units of ^l^y) and the wavelengths (shown in units of /.tm) are corrected to the 
rest-frame at the redshift noted. 



The IMF has a mass range of 0.1-125 M© and assumes 
solar metallicity. Higher metallicities are often favored 
when fitting the SEDs of Ultraluminous In frared Galax- 
ies (U LIRGS) in the local Universe (e.g. iFarrah et al.l 
VMm . but we found no appreciable difference in the es- 
timates of Lx/M for a range of likely values. 

We can combine the predicted i^-band luminosities as 
a function of time from Starburst99 with the inte- 
grated stellar masses at that time to derive the effective 
light-to-mass ratio for the stellar population: 

Li^/M = IQ-''-^'^^'''-^-^'^ /^{AUt)), (1) 

where 3.3 is the absolute magnitude of the Sun in the K- 
band l)Coxll2000|) . We find that the Hght-to -mass ratio 
for the CSF case is reasonably described by a simple 
power law, L/^/M= 103(T/Myr)^° ''8 for t > 100 Myr. 
The behavior of /M for a burst model is not as easy 
to parameterize over the likely range of ages of SMGs, 



so instead we use tabulated values. In solar units, we 
find that Lif/M is [2.5, 2.2, 13.3, 7.6, 4.2, 2.7, 1.9] at 
10, 20, 50, 100, 200, 500, 1000 Myrs respectively. We 
caution the reader that at ages of < 50 Myrs changes 
to the adopted IMF or metallicity can have significant 
effects on the Ijk /M and so the predictions for extremely 
young bursts are much more uncertain. 

In principle, we can constrain the ages of the dominant 
stellar population in each SMG from our SED model- 
ing. However, the correlated nature of the fits and the 
fact that we are likely sampling different mixes of stellar 
populations at different wavelengths due to differential 
reddening means that the results are highly uncertain 
for individual galaxies. To demonstrate this point, we 
note that all of the galaxies can be reasonably fit by 
either model star formation history, but with very differ- 
ent ages. This is a direct consequence of the fact that 
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an older galaxy with a constant star-formation history 
has a similar SED to that of a y o unger burst model, a 
point also noted in lShaolev et alJ l)2001ll2005j) . For this 
reason we prefer to use only ensemble averages for the 
derived parameters. If we restrict our hyper-z fits to 
the burst model only, we estimate ensemble averages of 
250 ± 190 Myr, Ay = 1.7 ± 0.2 and = -26.3 ± 0.3. 
Alternatively, restricting our fits to just the continuous 
star-forming model yields 2000±500Myr, Ay = 1.7±0.3 
and Mk = -26.4 ±0.3. 

For the average age from our fits of the burst model, we 
estimate a typical Lif/M=3.7, while the average for the 
continuous star-forming model is Li<-/M=2.6. Hence, the 
average stellar mass from the samples, estimated using 
Mk and L^/M, are essentially identical, log(Mburst) = 
11.3±0.4andlog(McsF) = 11.5±0.4. Based on these re- 
sults, we choose to simply use L^f /M=3.2 for all galaxies. 
Given the degeneracy between the SED fits, we use the 
average from the burst and CSF fit for each galaxy, 
and half the difference between them as a measure of 
the uncertainty. We summarize the results of our stellar 
mass estimates from the SED fitting in Table |21 The 
mean stellar mass is log(Mj,) = 11.4 ±0.4 for the sample 
of SMGs used here. 

3.2. The M*-Mbh relation 

Fig. 13 shows the rest-frame V-band and if-band lu- 
minosities fro m our SED fits a gainst the X-ray lumi- 
nosities from lAlexander "eTall (I2005bl^ . There IS an 
apparent trend between if-band and X-ray luminos- 
ity, which is reasonably fit by a power-law of the form 
hK (xL'^^''=^°-°^. To test the strength of the correla- 
tion, we applied a Spearman rank test. The correlation 
coefhcient is 0.68, with a 2.4tT deviation from the null- 
hypothesis that Li<- and Lx are not correlated. Remov- 
ing the three lower redshift sources in Fig.|2]does weaken 
the apparent correlation: the power law index becomes 
0.79 ± 0.32, while the Spearman test estimates a 1.8cr 
deviation from the null-hypothesis. While the modest 
dynamic range of our sample makes evidence for a corre- 
lation difhcult to find, as we discuss below, it does yield 
strong constraints on the normalization of the relation- 
ship. 

For both the full sample and with the three lower red- 
shift sources excluded, the trends with i^-band mag- 
nitude are three times tighter than that seen between 
the reddening-corrected rest-frame y-band magnitude 
and X-ray luminosities, yielding an RMS of 0.6 versus 
2.0 mags respectively. This difference most likely re- 
flects the increased sensitivity to reddening corrections 
in the rest-frame T^-band, in addition to the fact that 
near-infrared luminosities better probe the evolved stel- 
lar population. Since the uncertainty in the stellar mass 
estimates scale as 0.4Ai„ag dex, these results underline 
the advantages of the rest-frame K-haivA luminosities to 
estimate the stellar masses of these galaxies. 

If we assume the X-ray emission arises from accretion 
onto SMBHs, a s indica ted by the spectral analysis of 
lAlexander et al.ll)2005bD . we can interpret the correlation 
between Lk and Lx as one between the stellar compo- 
nent of the galaxies and their SMBHs. Such a correlation 
would also be expected if the if-band light were due to 
hot dust powered by the AGN. However, with only a few 
exceptions, the SEDs seem to be dominated by stellar 
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Fig. 2. — We show the correlation between the rest-frame V- 
and K-haiid luminosities and the rest-frame absorption-corrected 
2-10 keV luminosities of the 13 high-redshift, SMGs in our sam- 
ple. Open symbols represent the optical luminosities uncorrected 
for reddening. The best fit value of Ay from HYPER-z were used 
to correct these values for reddening, and we show them as the 
solid points. Circles represent the 10 z > 1.8 subsample, while 
triangles are the three lower redshift sources. The iC-band panel 
shows a clear trend for more luminous galaxies to have higher X- 
ray luminosities, as expected if they host more massive black holes. 
A similar trend may be seen in the V-band panel, however, this 
is masked by the increased sensitivity to the reddening correction 
and the corresponding uncertainties. In order to keep the scales the 
same between the two panels, we excluded J123712. 05+621212. 3 
from the V-band plot since it falls considerably outside the range 
(with My = -13.5). 



light, as evident from a well-defined stellar bump peak- 
ing near 1.6 /im(Fig. For the remainder of the paper 
we will assume that the if-band is tracing the stellar 
mass, but the reader should bear the caveat in mind. 

In Fig. O we plot the distribution of M^, - the 
stellar mass from the SED fits assuming a fixed 
Lx/M, and MeH.Edd - the SMBH masses derived 
from the absorption-co rrected X-ray luminosities by 
lAlexander et alJ l|20fl5a|) assuming Eddington accretion 
and a 6l£j% bolometric correction. As in Fig. [SJ we see 
a reasonable correlation between these two quantities. 

The limited dynamic range of our sample and the un- 
certainties in the individual measurements means we can- 
not provide a reliable constraint on the power-law index 
for this correlation, Mbh ocM^. The data are consis- 
tent with 7 ~0.3~1 with a very large uncertainty de- 
pending upon the exact choice of sample used. Nev- 
ertheless, while the slope of the correlation is uncer- 
tain, the normalization is well-constrained at the mean 
M* of our sample. At log(M*) = 11.4 ± 0.4, we find 
that the average Eddington-limited black-hole mass is 
log(MBH,Edd) — 6.86 ± 0.40. If we exclude the three 
SMGs at z < 1.8, the mean stellar and black hole 
masses are slightly larger: log(M^) = 11.5 ± 0.32 and 
log(MBH,Edd ) = 7.00 ±0.27 

As a comparison we turn to a sample of nearby 
(< 250 Mpc) galaxies with reliable black-hole masses es- 
timated from high resolution HST spectroscopy, and stel- 
lar masses derived via near-infrared observations of the 
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TABLE 3 

Derived optical and X-ray parameters for a sample of SMGs in GOODS-N 





-a 

Z 








lOg^iViBH j 


Fit quality 


J123636.75-I-621156.1 


0.557 


-24.10 ±0.06 


10.46 ± 0.02 


42.7 






J123721. 87-1-621035. 3 


0.979 


-26.14 ±0.03 


11.27 ±0.01 


43.7 






J123629.13-I-621045.8 


1.013 


-26.32 ±0.05 


11.34 ±0.02 


43.2 


6 3+°-^ 
"•■^-0.3 




J123555.14-I-620901.7 


1.875 


-27.28 ±0.37 


11.73 ± 0.15 


44.4 






J123711. 98-1-621325. 7 


1.996 


-25.48 ±0.40 


11.01 ± 0.16 


43.6 






J123635.59-l-621424.lf 


2.005 


-27.98 ±0.49 


12.01 ± 0.20 


44.0 


7 1+0.5 
'■^-0.4 


Near-IR Excess 


J123549.44-I-621536.8'' 


2.203 


-27.44 ±0.32 


11.79 ± 0.13 


44.0 


7 1+0.5 
'■^-0.5 


Near-IR Excess 


J123606.72-I-621550.7 


2.416 


-27.57 ±0.30 


11.85 ± 0.12 


43.8 


7Q+0.5 
'■'^-0.4 


Near-IR Excess 


J123622.65-l-621629.7f 


2.466 


-26.34 ±0.58 


11.35 ± 0.23 


44.0 


7 9+0.5 
'■^-0.5 




J123707.21-|-621408.lf 


2.484 


-26.40 ±0.19 


11.37 ±0.08 


43.8 






J123606.85-l-621021.4f 


2.505 


-27.00 ±0.37 


11.61 ± 0.15 


43.7 


6.8tE!:^ 


UV Excess 


J123616.15-I-621513.7 


2.578 


-26.60 ±0.12 


11.46 ± 0.05 


43.7 


6.8l^:^ 




J123712.05-I-621212.3 


2.914 


-25.82 ±0.15 


11.15 ± 0.06 


43.4 







Note. — 

^Redshift determined from UV spectroscopy IChapman et al .1200511 . 

^Absolute 7i"-band magnitude derived from HYPER-Z (uncorrected for reddening). 

^Stellar mass in units of Mq derived assuming Lk /M— 3.2. 

*^Rest-framc 0.5-8.0kcV luminosity (in units of ergs~^) and SMBH mass in units of Mq from [Alexander et alj l2005aD assuming a bolometric 
correction of 6%, with an uncertainty reflecting the possible range in this correction. 
^SED fit was deemed good unless otherwise notes. See §3.1. 

^These five SMGs also have rest-frame optical spectroscopy from lSwinbank et alj 1200411 




IQIO IQll 101=^ 



Stellar masses (M./Mg) 

Fig. 3. — The relation between stellar and SMBH masses for our 
sample of 13 high-redshift SMGs. Black hole masses are de rived 
using the assumption of Eddington-limited accretion (see i|4.2l . 
Stellar masses are derived by using the rest-frame ii'-band absolute 
magnitude derived from the fits to the SEDs in Fig.^and assuming 
Lx/M=3.2 for all galaxies. The upper dashed lin e denotes the best 
fit to l ocal galaxies (shown as stars) as derived bv lMarconi fc HunH 
<2003l) , and the grey area denotes the region where 68% of the local 
points are contained. The lower dashed lines denote tracks with 
black hole masses a tenth and hundredth of the local relation. The 
three SMGs shown as open symbols are those lying at redshifts of 
2 < 1.8. 



host galaxies (|Marconi fc Hunt! I2fl0,'^ . Although their 
full sample includes 37 galaxies, we concentrate on the 
27 "Group 1" systems judged b y those authors to be the 
most secure. Using these data, iMarconi fc Hunt) l)2003fl 



showed that: 

log(MBH) = (8.28 ± 0M)+ (2) 
(0.96 ±0.07) log(MBuige/10i°-9). 

We note that this relation predicts black hole 
ma sses roughly 50% higher than those presented in 
iHaring fc R ix ( 2004) , who used Jeans equation modeling 
rather than virial estimators to predict the bulge mass. 

Using these relations, our data suggests that the black 
hole masses we infer for SMGs (either in the full sample 
or for the ten SMGs at z > 1.8) are ~ 50 times smaller 
than for similarly massive local spheroids. 

4. DISCUSSION 

With the assumptions about the Lk /M and Eddington 
fraction t] we have adopted, our results suggest that the 
black holes in SMGs require significant growth in order 
to approach the present-day relation. To understand the 
significance of this conclusion, however, we need to better 
characterize any biases in the mass estimates for both of 
these components. Note that errors on individual sources 
are large in principle, since for both the black-hole and 
stellar masses we have assumed ensemble values (i.e. fixed 
77 for the black-holes, and a single Lk /M value for all 
stellar mass estimates). It is quite likely that source-to- 
source variations will affect the scatter of the data, but 
until more accurate methods are available for the mass 
estimates, we concentrate instead on understanding any 
systematic bias which could effect the sample as a whole. 

Altogether, the ratio between between our black hole 
masses and those in the nearby Universe scales as: 
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Age/Myr 
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K-band Light-to-Mass Ratio (l^/M) 

Fig. 4. — This figure graphically illustrates the available pa- 
rameter space required to place the SMBH and stellar masses of 
the high-redshift SMG sample onto the local M»-Mbh relation. 
The solid lines represent combinations of /M and ?j that would 
bring our mass estimates in line with the local relation at a given 
assumed value for the X-ray bolometric correc tion. They span the 
estimate for BCx determined by Elvis ct al] 119941) . We plot our 
best estimate for the values of L^/M and rj as the solid point. 
As an example to help illustrate the use of this figure: assume a 
6% X-ray to bolometric luminosity correction and an Eddington 
efficiency of 0.1, the plot indicates that our stellar ages need to be 
~ 50 Myr in order to make the data agree with the local relation- 
ship. The plot does not account for a possible missing fraction of 
X-ray fiux due to heavy obscuration in the SMGs. This would sys- 
tematically push the black-hole mass es timat es u p, and the three 
lines c loser to our observed point. See H4.2I and lAlexander et afl 
12005hi) for more details. 



where we have simplified the expression by assuming the 
power-law relation between stellar and black hole has an 
exponent of 1.0. In Fig.0|we plot tracks of different con- 
ditions where the above expression is unity. This demon- 
strates how difficult it is for any one factor to explain the 
difference between our sample and the local relation, and 
suggests that either the discrepancy is real, or that bi- 
ases in the parameters must conspire if, in fact, SMGs 
have a black- hole/stellar mass ratio comparable to local 
systems. We now discuss these potential biases in more 
detail. 

4.1. Factors influencing stellar mass estimates 

We begin by pointing out that the stellar masses we 
derived are on average five times higher than those re- 
po rted from modelin g the UV/optical SEDs of SMGs 
in iSmail et al.l l)2004fl . This contrasts with the reason- 
able agreement in the estimated stellar masses derived 
for UV-selected galaxies at similar redshifts to the SMGs, 
by Shapley et al. (2005) from both rest-frame UV/optical 
and rest-frame UV/optical/near-infrared photometry. 
We attribute the d ifference in our masses to those from 
ISmail et alJ 1)20041) to the absence of mid-infrared infor- 
mation to constrain the SED modeling in the earlier 
work, which is then much more sensitive to the redden- 
ing corrections derived from the SED fits to these very 
dusty galaxies and the possibilities of strong differential 
reddening within their young stellar populations. 



In our modeling of the stellar mass of the systems we 
adopted a Lk/M ratio of 3.2, characteristic of a burst 
with an age of ~ 250Myrs. If we require that the SMG 
and the present-day M^^-Mbh relations agree, then from 
Eq.Elwe see that the stellar masses must drop by a factor 
of ~ 50. This corresponds to L^/M^ 250, which is not 
possible with the Miller-Scalo IMF. For a Salpeter IMF, 
such a high Lk/M ratio implies ages for the SMGs of 
< lOMyrs. We discuss the implications for varying the 
IMF later in this section. 

However, we note that there are fundamental incon- 
sistencies in our calculations if we rigidly interpret the 
star formation histories we used. The SMGs in our sam- 
ple have been selected on the basis of their intense far- 
infrared luminosities, which implies high star formation 
rates. The fast decline in the far-infr ared emission fro m a 
starburst after star formation ceases l|Kennicuttll998)) in- 
dicates that star formation is currently on-going in these 
systems or has declined significantly in only the last few 
lO's of Myrs. The submm detection of these sources is 
therefore formally inconsistent with the star formation 
history in our burst model - where there has been no 
star formation within the galaxy in the past 250 Myrs. 
Although the alternative CSF models include on-going 
star formation, the much higher ages derived from these 
fits are equally inconsistent with a constant rate of star 
formation at the level indicated by the observed far- 
infrare d emission from the SMGs, ~ 1-2 x 10'^ M© yr~^ 
(jKenn icutt 199|). To form an Mk ^ -26.4 galaxy with 
a constant SFR at this level would take just ^ 100 Myrs, 
not the 2 Gyrs derived from the fits. 

There are two solutions to solve these inconsistencies 
within our models: 1) allow a more complex star for- 
mation history, most likely involving a recent burst of 
activity; 2) apply age-dependent dust extinction to the 
stellar populations, with the youngest stars being most 
extincted. As both solutions involve additional free- 
parameters to what is already a weakly-constrained prob- 
lem, we will only qualitatively discuss their impact. 

Bursts of star formation are a natural consequence of 
interaction-induce d act ivity due to the cyclic behavior of 
bound orbits l|Mihos fc H crnauist 1994). They lead to a 
dispersion in the properties of SMGs dependent on the 
number of previous bursts they have experienced. The 
main consequence in terms of our analysis is the realiza- 
tion that the present activity may not have produced the 
total stellar population we observe. Mixing in a fraction 
of older stars will systematically decrease the effective 
Lk /M, potentially increasing the masses of the galaxies. 
This would therefore push the SMGs even further away 
from the z = M^-Mbh relation. 

Age-dependent or selective ext i nction of the stellar 
populations (e.g. iPoggianti et al.l 1200 ID within SMGs 
can preferentially extinguish the UV emission from the 
youngest stars. This results in SEDs with much redder 
spectral indices between the UV and near-infrared than 
can be achieved with a simple foreground dust screen 
acting upon the whole stellar population. When the dif- 
fusion timescale for stars to escape their obscured birth- 
places is around a few 100 's Myrs, this process can also 
help reproduce the strength of the spectral breaks at 
^ 4000A, even in galaxies with significant on-going star 
formation. As a result, selective extinction would pro- 
vide the opportunity to fit the SED of our SMG sam- 
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pie with relatively young CSF models, ~ lOOMyrs, con- 
sistent with the ages needed to produce their observed 
near-infrared luminosities given the star formation rates 
inferred from their far-infrared emission. 

Attempts to estimate typical ages for SMGs include the 
spectral dating of the few example s with high-quality UV 
spectroscopy. Such an analysis bv lSmail et al.l l)2003(l in- 
dicates an age for a UV-bright SMG of ^ 10 Myrs, similar 
to the minimum age we require. However, as lSmail et al.l 
l)2003(l state the SMG they model may be a rare sys- 
tem and the small proportion of SMGs showing simi- 
lar signs of young stellar populations in their spectra 
suggests a more typical age for an SMG would be 50- 
100 Myrs. Another estimate of the likely age of SMGs 
can be derived via gas-depletion timescales using the 
measured mass of the gas reservoir and the estimated 
star formation rates suggest likely lifetimes of Z 30 Myrs 
ijGreve et al.l 120051) . where we have assumed an AGN 
contribution to the far-i nfrared luminosities in line with 
ijAlexander et al.ll2005btl . Finally, taking our estimated 
stellar masses and the median star formation rates de- 
rived from the far-infrared luminosities of the SMGs in 
our sample using our adopted IMF, we determine life- 
times of ~ 100 Myrs to produce the stellar populations 
we see. Balancing these estimates with those from the 
SED fitting, we conclude that the typical ages of SMGs 
are likely to be Z 50 Myrs. Hence the stellar masses are 
inconsistent with what would be predicted from the lo- 
cal M^^-Mbh relation assuming our black-hole masses are 
accurate. 

Are there other ways to change the effective Lk/M 
without invoking such young ages? One way to reduce 
it is to have a significant contribution to the rest-frame 
ivT-band light from the AGN. However the bulk of our 
sample, as evident in Fig.Q] are well-fit by a stellar com- 
ponent only. In the few cases where AGN cause the SED 
fits to diverge at the longer wavelengths, the estimated 
luminosity we derive at 2.2/Ltm is still dominated by the 
stars. We used the 8 fim excess (the ratio of the mea- 
sured 8 /im flux over that derived from the best fit model 
template) to quantify possible AGN contamination. We 
find that the z > 1.8 sources (excluding the 3 marked as 
having an obvious near-IR excess) have an average excess 
of 1.0, suggesting little or no contribution from an AGN. 
The remaining three sources however, have an average 
excess of 2.3, which means a nearly equal contribution to 
the rest-frame X-band from AGN as stars. Scaling down 
the stellar masses of these three galaxies to account for 
this results in only a ~ 0.1 reduction in the mean stellar 
mass of the sample, which is small considering the other 
factors we discuss. 

In addition to contribution of light from the AGN, 
we note that our analysis includes light from all of 
the stars, and not simply the bulge component used in 
deriving the local relation. While the high-resolution 
ffST/NICMOS imaging verifies that the rest-frame V- 
band light in SMGs are reasonably fit by a r^/^ sur- 
face brightness profile l|Borvs et al.ll2006() . some systems 
do show a significant disk component. Whether these 
stars are still in a disk component at the present-day 
or whether they would have been driven into the bulge 
through, for example, bar-instabilities, is impossible to 
verify, and so we flag this to the reader as a fundamental 



uncertainty of our analysis. Nevertheless, we note that 
any correction to the stellar masses due to a disk com- 
ponent is unlikely to change the inferred spheroid mass 
by more than a factor of 2-3 x for a small fraction of the 
SMGs. 

The I MF also has a signific ant effect on stellar mass es- 
timates. iBaugh et al.l 1)20051) argue that an IMF biased 
to more massive stars is required to reproduce the prop- 
erties of SMGs in current semi-analytic galaxy formation 
models, and invoke a power-law IMF with an index flat- 
ter than the canonical Salpeter ( Salpctcr 1955) slope of 
a ~ 2.35. We used Starburst99 to study the effect of 
this top-heavy IMF on our mass estimates, but find that 
the stellar mass estimates change by at most a factor 
of ~ 3 for the range of likely ages for SMGs. An ad- 
ditional factor of three in the light-to-mass ratio can be 
picked up by truncating the IMF at ~ 1 Mq , and there 
ar e different IMF p rescriptions (such as the one proposed 
bv lMaraston et all (2005) that has an increased K-hand 
contribution from thermally pulsating AGB stars) that 
can also increase the effective light-to-mass ratio. Al- 
together it is possible that these effects could conspire 
to decrease our stellar masses by up to a factor of ^ 10. 
However, we conclude that any bias present in the stellar 
mass is insufficient by itself to account for the difference 
between our results and the local relation. 

4.1.1. Tests of the stellar masses 

Are the high stellar masses we derive in conflict with 
other estimates of the masses of SMGs? 

Dynamical masses for SMGs from CO spectroscopy 
with millimeter interferometers indicate a median 
gas-l-stellar masses of log(Mdyn) H-l ± 0.3 within the 
central ^ lOk pc of typical SMGs with 850-^m fluxes of 
a few mJys fFrave r et al]IT998l IT999t iNeri et al.ll200l 
Greve et ciL 2Q05J. In addition, flve of the galaxies in 
our sample have dynamical mass estimates from Ha 
spectroscopy (jSwinbank et alJl200'H) . which all give con- 
straints on the masses of log(Mdyn) ~ 11.0 -11.3, again 
on ~ lOkpc scales. In a recent analysis, iGreve et alJ 
f2005) use the integrated CO line luminosities for 12 
SMGs to estimate an average contribution from gas of 
log(Mgas) — 10.5 ± 0.2 within ^ lOkpc for these systems. 
Thus taking into account the contribution from gas, the 
dynamical estimates are consistent with the presence of 
log(M*) 11 of stars within the central regions of the 
SMGs. 

As highlighted by ISwinbank et all l|2004j) the Ha dy- 
namical mass estimates derived from the Ha line widths 
for SMGs are flve times higher than similar constraints on 
UV-s elected galaxies at comparable redshifts (Erb et aQ 
I200I . To see if this difference is also reflected in the stel- 
lar mass estimates for the two populations, we ran our 
HYPER-z analysis on the sample of 72 spectroscopically 
conflrmed, UV-selected, and IRAC-detected, galaxies at 
z 2.3 presented in Shaolev et al. (2005) These UV- 
select ed galaxies (lAdelberger et all I2004t iSteidel et alJ 
l2003() make an ideal comparison set, as they are more 
modestly star-forming galaxies at roughly the median 
redshift of our s ample. As with our analysis, the flts 
to the data from IShaplev et all l|2005f) were degenerate 
in SED template, age and reddening. Roughly one-third 
of the objects are best flt by a young (~ 50 Myr) burst 
while the remainder are better modeled by a constant 
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star-formation history with age of ~ 600 Myr. Again we 
find that the choice of model has httle impact on the 
derivation of the steUar mass, and we obtain a mean of 
log(M,) = 10.35 ± 0.34 for this sample. The estimated 
ages and mas ses are in go od agreement to the results 
which Shaplev ^et al.l l)200!Tt) derive using a more detailed 
analysis (log(MH.) = 10.32 ± 0.51). This confirms that 
our approach to stellar mass estimation is robust and 
also suggests that the SMGs in our sample are on av- 
erage ~ 10 times more massive than the UV-selected 
population at z ~ 2.2. 

Overall, we conclude that our stellar masses are con- 
sistent with dynamical limits on the masses of the SMG 
population and support the suggestion that the central 
regions of these galaxies are baryon dominated. We note 
that the CO detections of SMGs indicate that there is 
sufficient gas present to continue to build the black hole, 
but at the same time the gas masses are small enough 
that they will not appreciably change the stellar masses 
even if the bulk of it is converted into stars. Only by 
bringing in more gas or stars into these systems (through 
cooling or mergers) can their stellar masses be signifi- 
cantly increased from what we measure. 

4.2. Factors influencing black hole mass estimates 

Is it possible that the sources only appear to lie below 
the l ocal relation because the black-hole mass estimates 
from [Alexander et al.l l|2005a(l are too low? Certainly, 
in the light of the discussion above it appears that the 
black hole mass estimates are probably the more impor- 
tant source of uncertainty in o ur analysis. In e stimating 
the black- hole masses. [Alexander et al.l l|2005aj) assumed 
Eddington-limited accretion and a 6^4^% bolometric cor- 
rection to convert the observed absorption-corrected Lx 
to MBH,Edd- Since the Eddington limit appears to pro- 
vide a reasonabl e upper bou nd to the accretion rate of 
black holes (e.g.. lMcLure fc D unlop 2004), this approach 
essentially provides a lower limit on the black-hole masses 
and sub-Eddington accretion would imply larger black 
holes. However, an average black-hole mass ~ 50 times 
higher is needed to reconcile our data at the character- 
istic stellar mass of the SMGs, ~ IO^Mq. This would 
require rj ~ 0.02, which is low given the gas-rich environ- 
ments in which the SMBHs are likely to reside. Although 
we do not know the Eddington rate of these galaxies, the 
large fraction of SMGs that host AGN activity suggests 
that the accretion is like ly to be reasonably efhcient (see 
I Alexander et al.l lj2005aj) '). There is tangential support 
for this claim from studies that have estimated black- 
hole masses for other active high-redshift populations 
independe ntly of the Eddington r ate, findin g rj ^0.1- 
1.7 (e.g. E larconi ct al. 2004; McLure fc Dunlorii [260I 

iBarger et al 2000.') . At the smallest accretion efficiencies 

found in these studies, the black-hole masses would still 
be 5 times too small for our estimated stellar masses. 

There is further support for Eddington-limited or 
near Eddington-limited accretion from the similarity 
between the optical/near-infrared spectra of some SMGs 
with thos e of local Narrow Line Seyfert 1 (NLSl) 
galaxies (llvison et al.l 1199^ iVernet fc Cimattil 1200 It 
Ledlow et al.l I2002t ISmail et alJ 1200 3: Daws on et all 
2003: "Swinbank et al.' 12004 I2005|) : see §4.3 of 
Alexander et ah. l,2005b|) . If the spectral similari- 
ties of the AGN signatures in these two populations 



are suggestive of comparable physical conditions in the 
accretion disk around the SMBH, then the approx- 
imately Eddington-limited rates estimated for local 
NLSls may be appropria t e for SMGs (e.g. iBorosonI 
l2?ffll ICoUin fc Kawaguc hi' "2004^. Moreover, NLSls 
appear to have smaller black holes than those predicted 
from the l ocal M^,-Mrh rel at ion; U sing r everberation 
mappm g, iGrupe fc Mathud (|2004') and iBotte et alJ 
(j2004|) show that that the SMBHs in NLSls are an 
order of magnitude smaller than predicted for galaxies 
having a stellar mass of lO^^M©. Thus physical models 
for the accretion disks and SMBHs in NLSls may be 
very relevant for understanding the growth of SMBHs 
in massive, young galaxies at high redshifts. We do 
caution the reader that the broad-line region in SMGs 
may be missed due to heavy extinction in these dusty 
galaxies, and hence the spectra may be dominated by 
a less obscured narrow-line region. High signal-to-noise 
spectra would be useful to search for faint, broad wings 
that may be otherwise missed. 

There is also theore tical supp ort for ri ^ 1 efficie ncies: 
iGranato et"al] l)2004|) and Kaw akatu et all l)2003|) pre- 
dict that this condition will be met near the peak epoch 
of star formation, where these objects have a stronger 
chance of being detected at submm wavelengths. In- 
terestingly, these models also predict black hole masses 
« 1-2 orders of magnitude smaller than the local rela- 
tion during the obscured, star-bursting phase, similar to 
what our data suggest. Hence, while the assumption of 
Eddington limited accretion is a large one, it is not a un- 
reasonable choice. We reiterate that even an efficiency 
10 times smaller would still result in the black-holes in 
SMGs being smaller than would be predicted by the local 
relation. 

Another uncertainty in the SMBH estimates comes 
from the bolometric corre ction applied to the X-ray lumi- 
nosities. [Alexander et alJ l|200 5a') assumed that 6^4^% of 
the energy out put due to ac cretion was emitted in the X- 
ray wavelaand ijElvis et alJ lT994). but if this was lower, 
the inferred mass of the SMBH would rise. However, 
they also pointed out that this correction is similar to 
the estimated bolometric corrections for obscured and 
unobscured A GN of similar luminosity to the AGNs in 
the SMGs (see [Alexander et al.ll2005a() . 

A final potential source of uncertainty in the SMBH 
masses comes from the possibility that a substantial 
fraction of the X-ray emission from the SMBH is still 
obscure d at rest-frame ene rgies of ~ 20 keV. For in- 
stance, iWorslev et all l)2005ft find that the Chandra 2- 
Ms exposure of the CDF-N fails to recover approxi- 
mately 40% of the hard X-ray background. This may in 
part b e hidden within highly-obscured 'Type-2' quasars 
(jMartinez-Sansigre et al. 2005), but dusty SMGs also 
harbor obscured AGN (and could in some cases be Type- 
2 quasars). Indeed, not all SMGs are detected in the 
Chandra 2-Ms X-ray image and so potentially some frac- 
tion of the X-ray emission from the SMBH in this popula- 
ti on could be missed from the Chandra analysis. In §3.5 
of I Alexander et all l)2005bj) . it is pointed out that while 
obtaining tight constraints on the obscuring column den- 
sities for individual sources is difficult, the agreement be- 
tween their models and X-ray spectra of SMGs ranked 
by obscuration class is excellent, suggesting that most of 
the X-ray emission in these sources is accounted for. If 
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a large fraction of the X-ray emission in SMGs is still 
absorbed, we will need to wait for the next generation of 
X-ray observatories to find it. 

4.3. Implications for the evolutionary history of SMGs 

To summarize the previous two subsections, it is pos- 
sible to find combinations of the parameters used to es- 
timate the stellar and black hole masses of SMGs which 
would reconcile them with the local M^^-Mbh relation. 
For example, if we assume a 6% bolometric correction, 
then we require a combination of a typical SMG age of 
just ~ 50 Myrs and an Eddington efficiency of 77 ~ 0.1 to 
place the high-redshift data on the local relation (Fig. ^ . 
This combination is plausible, and if it is the case then 
the implication is that the stellar and black hole mass 
grow simultaneously, with a mass ratio that is similar at 
high redshift as well as low. However, the need for both 
low Eddington accretion and a very short lifetime for the 
SMGs is uncomfortable and hence we should also explore 
the possibility that the black holes in SMGs are, in fact 
less massive than expected from their stellar masses and 
the low redshift M^^-Mbh relation. The following discus- 
sion assumes that this is ind eed the case. 

The theoretical models of iGranato et al.l |2004) sug- 
gest that SMBHs may grow exponentially, taking only 
a few IP's of Myr to mature ( Kawakatu et al. 200^ 
IGranato et al]l200l . This is a short enough time-scale 
that SMGs with more evolved SMBHs, which lie on the 
local relationship, should fall within the redshift range 
probed by our sample. That we do not see such sys- 
tems implies that as the black hole rapidly builds up its 
mass, it either shuts down the star-formation that pow- 
ers the far-infrared luminosities (hence they drop out of 
the submillimeter samples) or they become so obscured 
that they are no longer detected in the X-ray waveband. 

The former argument supports a scenario that links 
the QSO and SMG populations. As the black hole un- 
dergoes a rapid buildup of mass, the ensuing strong ra- 
diation field strips away the obscuring gas and dust, 
terminating star formation and eventually allowing the 
AGN to be visible for a short time as a luminous opti- 
cal qu asar (Sanders et al. 1988; Silk & Recs 1998; FabianI 
119991 loTMatteo et ah .2005; .Alexander et a l. 2005^ 
The close similarity between the redshift distr ibution of 
radio-identified SMGs and QSOs ijChapman et al.ll20d5^ 
also lends support to this hypothesis, as does the relative 
space densities of SMGs an d QSOs given our estimates 
of their respective lifetimes (|Smail et alJl2003) . In addi- 
tion, recent high resolu tion adaptive optics imaging by 
IKuhlbrodt et alJ l)2005j) has shown that the black hole 
masses and host galaxy luminosities of z ^ 2 quasars 
are consistent with the local M*-Mbh relation, suggest- 
ing that QSOs could be the end-product of the rapid 
growth of SMB Hs seen soon after an S MG-like phase. 
Finally, work bv iPaEfe et all l|20?)ll I200I using SCUBA 
to measure the rest-frame far-infrared luminosities of un- 
obscured and moderately obscured QSOs at z ~ 2 ver- 
ify that the former class represent a phase of significant 
SMBH growth but not bulge growth, while the latter 
may be a transition phase between the star-formation 
dominated SMGs and AGN-dominated quasars. Thus, if 
SMGs are the progenitors of quasars, they must lie below 
the locally determined M^,-Mbh relation. 

In conclusion, there is good qualitative and quantita- 



tive agreement for an evolutionary sequence which con- 
nects moderate mass SMBHs in far-infrared luminous, 
star-forming galaxies at z ~ 2-3 (when most of the 
spheroid is produced), with subsequent growth of the 
SMBHs during a short phase as a luminous quasar to 
result in massive spheroids at the present-day which lie 
on the M^-Mbh relation. 

5. CONCLUSION 

Using deep Chandra, Spitzer, and ground-based op- 
tical and near-infrared imaging data, we have estimated 
the stellar and black hole masses for a sample of 13 SMGs 
with precise spectroscopic redshifts. Our main conclu- 
sions are: 

1. There is a tentative correlation between the rest- 
frame near-infrared and X-ray luminosity in this 
sample of SMGs. The X-ray emission likely arises 
from accretion of matter onto a SMBH, and so can 
be used to estimate its mass. We find a similar 
correlation between rest-frame X-ray and V^-band 
luminosities, but this exhibits more scatter, sug- 
gesting that it is more sensitive to dust and age 
variations in these young, obscured galaxies. 

2. The restframe UV-optical-near-infrared SEDs of 
most of our SMGs are dominated by stellar emis- 
sion and not light from the AGN, and are well- 
described by either a young burst model, or by an 
older galaxy with a more continuous star-formation 
history. In both cases, the derived light-to-mass ra- 
tios, and hence stellar mass estimates, are similar. 

3. The mean stellar mass for the sample is log(M^,) = 
11.4 ± 0.4, consistent with Ha and CO dynamical 
mass estimates for the SMG population. Hence the 
Spitzer observations have shown that these intense 
star-forming galaxies have a significant stellar pop- 
ulation in place at z ~ 2.2. 

4. We find that SMGs are typically 10 times more 
massive than typical UV-selected galaxies at simi- 
lar redshifts. 

5. We use the assumption of Eddington- limited ac- 
cretion to estimate the black hole masses for each 
galaxy in the sample. We find that for a given stel- 
lar mass, the black holes are ~ 50 times less massive 
than those for spheroids with similar stellar masses 
in the local Universe. 

6. In order to determine if the high-redshift SMGs 
lie of the local M^-Mbh relation, we explore the 
range of assumptions that are used to estimate 
both the stellar and black hole masses. We find 
that most reasonable combinations of parameters 
lead to lower black hole masses at a fixed stellar 
mass in SMGs at z ^ 2 compared to the low red- 
shift relation for spheroids. 

7. We suggest the NLSl galaxies may be a good tem- 
plate for understanding the growth of SMBHs at 
high redshift. They have spectral properties com- 
mon with many SMGs, and (locally) have black 
hole masses smaller than predicted for less active 
spheroids. 
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8. Our results arc consistent with a model where 
the SMBHs in SMGs subsequently undergo rapid 
growth to reach the local M^-Mbh relation. Such 
a phase would most naturally be explained if the 
SMGs evolved into quasars on their way to becom- 
ing luminous, passive ellipticals at the present-day. 
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